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ABSTRACT 

In this paper we first discuss the effects of oxidation of Si which estab- 
lished that vacancies (V) and Si self-interstitials (I) coexist in Si at high 
temperatures under thermal equilibrium and oxidizing conditions. Some essential 
points associated with Au diffusion in Si are then discussed. Analysis of Au 
diffusion results allowed a determination of the I-component and an estimate of 
the V-component of the Si self-diffusion coefficient. A discussion of theories 
on hign concentration P diffusion into Si is then presented. Although presently 
there still is no thoery that is completely satisfactory, significant progresses 
are recently made in treating some essential aspects of this subject. 


1. INTRODUCTION 


Atomic diffusion in Si can proceed via both direct and indirect mechanisms. 
Impurity atoms having no strong bonding interactions with Si atoms, and are 
hence located exclusively in interstices, can jump directly between the + inters- + 
tices. Species such as H, He (and presumably other noble gases), and Cu and Au 
(presumably the Cu and Au interstitials, respectively) are supposed to diffuse 
in this manner. Oxygen atoms, though possess strong bonding interactions with Si 
atoms, are also believed to diffuse directly by jumping between the bond-centered 
interstitial positions. By contrast, substitutional impurity and self-diffusion 
In Si need intrinsic point defects as diffusion vehicles he vacancy (V-) 

mechanism is known lo control ^elf-dif fusion in metals a. n Ge, and, as will 

be discussed, in Si below about 1000 C. Above 1000 C, however, the interstitial- 
cy (I-) mechanism plays a prominent role in Si self-dif f osion as well as in the 
diffusion of substitutional dopants P, B, A1 and Ga. 

The nature of the dominant point defect species in Si at high temperatures 
has been a long controversial subject, though it i<= clear that there can only be 
three possibilities: V only, 1 only, and I and V coexist. The controversy arises 
for two reasons: (i) Presumably because of the small point defect concentrations 
in Si, direct experimental methods have not been helpful. There is no reported 
absolute measurement (1) results, and there s only one very recently reported 
quenching result (2) which indicated the presence of V (in the opinion of the 
pr ent authors, however, th.s dees not imply that I do not exist); (ii) Theore- 
tical calculations, which are mere suitable for low temperat e Coses, have not 
yielded the needed ~5eV activation enthalpy and the -10k (k: Boltzmann constant) 
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activation entropy commonly observed for Si self— diffusion. There were several 
early discussions associated with some aspects of the possibility of I and V 
coexisting in Si: Seeger and Chik (3) considered that 1 and V may be contribu- 
ting equally to Si self-diffusion at 800°C; Prussin (4) considered that I and V 
should reach a local dynamical equilibrium under oxidizing conditions; Hu (5) 
considered that 3 should be diffusing in Si via an I -mechanism while As via a V- 
mechanism at high temperatures. In the last few years progesses in two areas 
helped to clarify our understanding of the nature of point defects and diffusion 
processes in Si considerably. Analyses of oxidation effects on stacking fault 
(SF) growth kinetics (6,7) and on dopant diffusion (8-10) showed beyond reasona- 
ble doubt that I and V coexist in Si at high temperatures under thermal equili- 
brium as well as oxidizing conditions. Analyses of Au diffusion in Si showed 
that I must exist (11,12) and that the features are totally consistent with the 
idea that I and V coexist (13) . Studies of Au diffusion into dislocation-free 
Si allowed to determine the I-component (14) and to estimate the V-eomponent 
(15) of the Si self-diffusion coefficient. In this paper the essential points 
associated with these two areas are first discussed. We then examine theories on 
high concentration P diffusion into Si. For various reasons, we believe that 
there is still no satisfactory theory. However, progresses have been made in 
treating some essential aspects of this subject. 


2. OXIDATION EFFECTS 


Oxidation of Si surfaces leads to the so called oxidation-enhanced and 
-retarded diffusion (OED and ORD) of the substitutional dopants and to the gene- 
ration of the oxidation-induced stacking faults (OSF) . These phenomena need to 
be considered together and are due to ^he fact that oxidation injects 1 into Si 
(5,16). It has been shown that a consistent interpretation of the data on OED/ 
ORD phenomena is not obtained if either I or V alone were assumed to be present 
under thermal equilibrium conditions (10). Therefore for our present purpose, 
we consider that I and V are both present in Si at nigh temperatures under ther- 
mal equilibrium conditions. This requires that we consider I and V have attained 
a local dynamical equilibrium condition (4,17), which was actually fulfilled for 
long time experiments (8-14,18,19). The condition is given by (17) 


C C 
I V 




( 1 ) 


where C^. and C^ r denote the I and V concentrations respectively and the supersc- 
ript eq denotes thermal equilibrium values. Eq . (1) is arrived at via the 
reaction 


T + V £ C 


( 2 ) 


(where 0 denotes the ideal lattice) under steady state conditions. Using some of 
the more reliable OSF size data (20-22), phenomenological but satisfactory ana- 
lyses of the OSF growth/shrinkage kinetics have been carried out (6,7). Define 
the I and V supersaturation ratios respectively as S^-C^/C^-l and 

S tr“C /cf^-l during the oxidation, we obtain (7) 

V V v 
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for the 0|j? wi|h a radius R . Here A is the^^rea^per atom in the fault 
(6.38x10 cm ) , ft the atomic volume (2x10 cm ) and y the stacking fault 
eaergy (0.026eV/atom) . The dimensionless quantity a f has a numerical value of 
=2 and contains all factors related to interaction potentials between I and the 
Frank partial dislocation binding the OSF. Emperically, for (100) Si wafers 
oxidized in dry C^, the OSF sizes are fitted satisfactorily by (6,7) 


Rgp = 1640exp(- 


2.5eV . 3/4 
kT ;t 


4.8xl0 9 

kT 


«xp(- 


5.02e V 

kT 


)t 


cm. 


(4) 


We use Eq. (1) to relate S T and S by S =-S T / (1+S T ) . A comparison of Eqs. (3) 
and (4) then yields (8,9), 

S° 2 = 6.6xl0 -9 exp(2.52eV/kT) t" 1/4 , (5) 

where the superscript o2 denotes the fact that this quantity is only suitable 
for (100) wafers oxidized in dry 0^* Eq. (5) is a most important piece of infor- 
mation obtained from OSF studies: it is used to obtain quantitative fittings of 

0ED/0RD data. S is an instantaneous value, the time averaged v; ] ae used in 
the actual OED/ORD data fitting is obtained by replacing the pre-exponential 
value in Eq. (5) by 8.8x10 and by noting that t then denotes the oxidation 
duration. 


In the presence of both I and V, the substitutional dopant diffusivity is 
given by 

D S = Dj + D®, (6) 

s s 

where and are respectively the I- and V-component of the dopant diffusi- 
vity. Under an oxidation which perturbs the thermal equilibrium I and V concen- 
trations, the dopant diffusivity changes to 

D® = D®(C ] ./Cj q ) + D®(C v /C® q ). (7) 

s s s 

In terms of the normalized diffusivity enhancement defined as A =D /D -1, the 
use of Eqs. (1) and (7) yields (8,10) ° K ° X 


A S = (2G +G S -1) / (1+S ) , (8) 

ox III I 

where G ^ is the fractional I-component of the dopant diffusity under thermal 
equilibrium conditions defined as G^=D /D . Eq, (8) applies to & an instantaneous 
diffusion enhancement. In experiments tne time averaged value A is measured, 
but Eq. (8) still holds to a good approximation by having useS together with 
A A plot of Eq. (8) is shown in Fig. 1 for three values. The value of A^ 
is X either positive or negative (0ED or 0RD) depending on the value of G^. and 
S . In Fig. 2 we show a fitting of available Sb 0RD data to Eq. (8). It is seen 
that the fitting is quite satisfactory on i quantitive basis with G^=0.02. This 
kind of good fitting indicates that the model of I and V coexisting and attained 
local dynamical equilibrium is correct, particularly in light of the fact that 
no other reasonably thought of model can be equally satisfactory (10) . The use 
of Eq . (8) allows to determine G^ fo- a dooant at a given temperature. Some such 
values have been given elsewhere (15 # 23,24*. To briefly summarize, it is found 
that the I -mechanism plays a prominenc rol *• in the diffusion of B, Ga, P and A1 
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at high temperatures, i.e., G >0.5 for these dopants at a temperature above 
1100°C. 


The above discussion applies for experiments with sufficiently long oxida- 
tion times. For short oxidation times, a transient exists for I and V to reach 
the local dynamical equilibrium. This is due to the presence of an energy barrier 
against I and V recombinations. The existence of such an energy barrier was sug- 
gestec by Wertheim (25) to explain the fact that the rate of A-center (E -0.16eV) 
production by high energy electron irradiation is strongly dependent upon the 
irradiation temperature (below room temperature). For high temperature cases, Hu 
(26) suggested that such a barrier of a few to several eV may exists, so that to 
a large extent I and V may behave independently . The experiments on ORD of Sb 
(18,19) shown in Fig. 2 demonstrated that for oxidation times longer than a few 
hours local dynamical equilibrium according to Eq. (1) is reached. By performing 
similar experiments for shorter times (5 to 60 min.) Antoniadis and Moskowitz 
(9) showed that Eq. (1) is not readily satisfied when oxidatio- first started. 

At 1100 C they found a small OED for a 5 min. oxidation which gives way to ORD 
for an oxidation longer than 10 min., and, the diffusion retardation approaches 
the value expected for local dynamical equilibrium at about 60 min. Their fin- 
dings may be rationalized as follows: during oxidation I are injected into Si by 
the SiO^-Si interface but it takes about 1 hr at 1100°C for I-V recombination to 
proceed to such an extent that local dynamical equilibrium is reached. As a con- 
sequence, at the oegining there is no V undersaturation so that the oxidation 
^ced I supersaturation enhances Sb diffusion via the small I-component of its 
ifusivity (G^-0.02). 

We use Waite's theory of diffusion controlled reaction (27) to obtain from 


T dlf 5 r - /4,, “ SDr iv 


(9) 


an estimate jf the time requited to establish the local dynamical equilib- 

rium in the Antggiadij-Mosl^owi t j experiment if no recombination barriergWere 
present. With D -10 cm sec ana the recombination radius r^^bxlO cm, 

Eq . (9) yields . ^.<0 .05sec which is about 10 shorter than observed. We there- 
fore conclude tha£ I-V recombination is controlled by the overcoming of a recom- 
bination barrier that exceeds the Gi^bs free energy of diffusion (of I and V) by 
AG. In this picture the factor of 10' by which the experimental x value differs 
from that estimated from Eq . (9) arises from the Boltzmann factor exp(AG/kT). 
Antoniadis and Moskowitz interpreted their observation in terms of an enthalpy 
barrier corresponding to AH-1.4eV assuming that in the expression 


AG - AH - TAS (10) 

the entropy contribution is negligibly small. By contrast, Gcisele et al. (28) 
proposed that the main part of AG originates from the term TAS where AS is nega- 
tive, i.e., from an entropy barrier, whigj^ Is due. a consideration of the large 
pre-exponential factor associated v* th D 

Irrespective of the detailed origin of the barrier, with AG-1.4eV, the bar- 
ri?* limited reaction time, t for establishing the local dynamical equilib- 

n between I and V is given 6y 
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(11) 


t, ~ r T _ r )exp(1.4eV/kT) . 

Dar lv 

Eq. (11) is plotted In Fig. 3 using D^«0.5iD C^+D^C®**) * ith the quantities 

and taken from Eq. (21) which will be discussed. It is seen that T ^ ar 

increases rapidly with a decrease of the temperature and attained values of 2, 

17 and 170 days at 1000, 950 and 900 C respectively Even if it is assumed that 
these results are over-estimated by a fa 'tor of 10, values arc still 5hrs, 

l*7days and l?days at 1000, 950 and yuO C respectively. One significance of this 
finding is that in the expression 


D %r - D I (C I /C r> 


+ D v <c v /C v ,) ' 


( 12 ) 


the quantities C and C^, arising from experimental perturbation of the respec- 
tive thermal equilibrium values, are not related by Eq . (1) iu in experiment 
performed in such low temperature range and having time durations much less than 
the appropriate t values. For device processing in the temperature range of 
850 to 950 C for a few hrs, should there be any reason that the point defect 
thermal equilibrium values are perturbed, the change only applies to the pertur- 
bed species independently, i.e., I and V need to be considered as not having 
attained the local equilibrium condition. 


The fact that I-V dynamical equilibrium is attained extremely slowly at 
lower temperatures offers the possihlity of measuring the quantities and D^ T 
independent of each other and independent of the I and V concentrations. There 
already exist some 0SF and 0ED/0RD experiments conducted with the intention of 
measuring at 1100 C or above (29-31). However, because of the temperatures 
and oxidation times involved, I-V local dynamical equilibrium should have been 
reached in these experiments and we have shown that such experiments yield an 
effective diffusivity given as (15,24,32) 


d i” ° »,<?-#>/<<?♦<?>. <u> 

eq eq 

instead of D . Depending upon the relative magnitude of C^. and , which are 
unknown quantities themselves, Eq . (13) may yield an approximate value for either 
Dj. or D . We have argued that the diffusivity values obtained from these experi- 
ments (29—31 should be interpreted as (24,32). However, for the lower tem- 
perature range, e.g., 850 to 950 C, the I-V local dynamical equilibrium w:' .1 
not be reached for an experiment involves less than a day’s time. For such ^ases, 
effects associated with dopant diffusion or stacking fault growth would be solely 
due to the point defect species whose thermal equilibrium concentration has been 
perturbed by the experimental condition. This means the other species does not 
contribute to the 'mental deviations observed if its thermal equilibrium 

concentration has nc jeen also directly perturbed by the experimental condition. 
We can now suggest an experiment to directly determine D^: Hign concentration P 
diffusion into Si proceeds efficiently in the temperature range noted above and 
the diffusion causes a quite high I supersaturation in the Si interior. Thus, 
experiments similar to those performed for OED/ORD and OSF (29-31) can be carried 
out at the low temperature range (850 lo 950 C) with surface P diffusion repla- 
cing oxidation to obtain . Similarly, D can be determined upon finding an 
interface reaction which injects V into Si. Presently we do not know with 
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certainty of such a reaction. 


3. DIFFUSION OF GOLD INTO DISLOCATION-FREE SILICON 


Gold atoms, which under thermal equilibrium conditions may incorporated 
in Si either as substitutional (Au^Au^-) or interstitial (Au^:Au ) atoms, may 
diffuse via either the kick-out mechanism (11) or the Frank-Turnbull mechanism 
(33) . The two mechanisms have in common that long range transport of Au atoms 
occur via migration of Au. atoms which may either jump from interstice to inter- 
stice to remain as Au. atoms or from interstice to lattice positions to become 
Au atoms. This is true because though the thermal equilibrium concentration of 
Au S is much larger than that of Au., the mobility of Au. is much much larger than 
thit of Au and hence long range transport of Au may be ignored. However, the 
two mechanisms differ in the ways Au^ and Au g atoms interchange. In the kick-out 
mechanj the interchange involves I according to (11) 

Au, J Au + I. (14) 

i s 

This mechanism creates an I supersaturation in the crystal which is balanced by I 
out-diffusion. This means if the mechanism is operative than I are involved in 
Si self-diffusion but it does not a priori mean that V are not contributing to 
Si self-dif fusion. In the Frank-Turnbull mechanism the interchange process invol- 
ves V according to (33) 

Au. + V i Au . (15) 

l s 


This mechanism creates a V undersaturation in the crystal interior which is 
balanced by V in-diffusion. If operative, it means that V are involved in Si 
self-diffusion but it does not mean I are not contributing to Si self-diffusion. 
We now assume that local dynamical equilibrium betweem I, V, Au . and Au is 
established. The necessary and sufficient condition for this assumption to be 
true is that the three reactions (2,14,15) have reached their equilibrium state. 
This requires that two out of the three reactions to be sufficiently fast. It is 
not required that the third reaction ba equally fast, since if equilibrium con- 
ditions are satisfied by two of the reactions, then the third is automatically 
under equilibrium conditions. Under this assumption the normalized Au concent- 
ration C=C ./ c et * to be measured in an experiment on Au dif f usion into disloca- 
tion-free Si" may be approximately described by 


1C _ _i (D 

at ax k eff ax' 


with the effective diffusion coefficient given by (13) 
D 


^ = < c “ 2+ v>i c r>vi q/c r- 


In Eq . (1/) contributions to the Au^-Au^ interchange process due 
out and Frank-Turnbull mechanisms are accounted for. In deriviug 


(16) 

(17) 

to both kick- 
Eq. (17) the 
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following approximations are used 



(C^ : Au^ concentration) 


(18a) 




C » 
s 


(c^) 1/2 . 


(18b) 

(13c) 


If only V and no I are present Eq. (17) yields a constant effective Au diffusity 
(11,33) s 



D C eq 
V V 



(19) 


The same result holds if is not zero but among the reactions (2,14,15) only 

(15) operates sufficiently fast. If only 1 and no V are present Eq. (17) yields 
a strong concentration dependent diffusivity 


D eff " < 20) 

eq 

The same result holds if is not zero but among the reactions (2,14,15) only 

(14) operates sufficiently fast. As observed by Stolwijk et al. (14), Au concen- 
tration profile after diffusion into dislocation-free Si at and above 800 C can 
not be describe! by the constant diffusity given by Eq. (19) but are satisfac- 
torily fitted by the use of Eq. (20). Fig. 4 shows a typical experimental profile 
due to Stolwijk et al. It can be seen that an erf c-func|ion type profile as expe- 
cted for D failed to fit the data while the use of D is quite satisfactory . 
This shows that I are contributing to Si self-diffusion But it does not mean that 
V are not contributing for the reasons that: (i) It can not be certain that this 
is not only due to the fact that of the reactions (2, 1^,15) only (14) operates 
sufficiently fast; and (ii) The strong dependence of D ^ on C reduces the 
effect^of D by this factor for most part of the experiment for which C<1 
and C <<1 nold. Indeed Moreheao et al. (13) found that while almost perfect fit- 
ting were obtained for short time diffusion profiles Laiall but observable devia- 
tions occured if Eq. (20) were used but a much better fitting is obtained by the 
use of Eq. (17). This not only shows that V are contributing to Si self-diffusion 
but also allowed to estimate that D C^D at 1000°C. The quantity has 

been determined by Stolwijk et al. (14) as 


D = 914exp(-4 .84eV/kT) cm^sec \ (21a) 

i j~ 

We have already shown that at T<900°C reaction ( 2 ) is expected to become ineffec- 
tive. At still lower temperatures, as will be discussed, reaction (14) is also 
expected to become ineffective. In that case only reaction (15) will be operating 
and an erf c-f unction type profile should show up. Willcox et f °u^ci such 

a case at 700°C. With their data we ootained that D^C^^8.8xl0 cm sec at 
700°C (15). This together with the data of Morehead et al. (13) yields fae rough 
~ jtimate 
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(21b; 


DyC^ = 0.6exp(-4.03eV/kT) cm^sec 

Lq. (21) is plotted in Fig. 5 in which a data point on D obtained by analy- 
sing the Ni precipitation behavior (15) of Kitagawa et al.(35) is also shown. 

An important point is, as rough estimates, the sum of the formation and migration 
entropies are now obtained from Eq. (21) as 3k for V and 10k for I. That is, 
wjjjie I are quite extended (3), V are fairly point-like. Experimentally observed 
D ^ values from relatively higher temperature measurements are usually of 
10 exp(-5eV/kT) while those from relatively lower temperature experiments are 
of lexp(-4eV/kT) (36-40). These values are quite compatible with Eq. (21). 

We now make an estimate to show that at very low temperatures the kick-out 
mechanism, reaction (14), can become ineffective, while the Frank-Turnbull 
mechanism, reaction (15), still operate-. It is important to note that if this 
can not be true then it can not be explained why Au diffuses into dislocation- 
free Si apparently via the kick-out mechanism at higher temperatures but the 
Frank-Turnbull mechanism at a very low temperature. Eqs. (19) and (20) show that 
Au^dif fusion^would appear as dominated by the kick-out mechanism in long as 
C Vr/Dv C v q >l Because for a large part of the experiment C >>1 holds, 

the Rick-out mechanism may still appear as dominating the diffusion process even 
for a case for which the frequency of the kick-out events is less than that of 
the Frank-Turnbull events. On the other hand, the frequency of the kick-out 
events should decrease more rapidly than that of the Frank-Tuenbull events with 
a decrease in temperature because the former J nvolves the creation of I which is 
charaterized by an activation enthalpy (the I formation enthalpy) while for the 
latter there is not yet a good reason for assuming that a high enough energy 
barrier exists to prevent an Au. atom drops in a V. In the following we show that 
under thermal equilibrium conditions a quantitative estimate can be made which 
indicates that the above discussed situation can be true. 

For the Frank-Turnbull mechanism, if an Au. atom and a V meet, there is a 
definete probability T<1 that reaction (15) will go forward to produce an Au 
atom. Under thermal equilibrium conditions, we obtain the rate of increase in C 
via this mechanism as 

(dC s /dt) FT = r(4TTr iV /n)D.C* a ^ q . (22) 

For the kick-out mechanism, reaction (14), the Au^ atom has a certain frequency 
v by which it produces Au and I . We do not know this frequency but if Au and 
I^meet there is a def inite S probabi lity 0<1 for reaction to produce Au^ and S under 
thermal equilibrium conditions this reaction balances that due to v^. Therefore 
we have , 

(dc /dt), = 0(4irr /i7)D T C® q C eq . (23) 

s ko Is I I 3 

From Wilcox et al. (34) * 3 have 

D , C*f q = 1.8xlCf 2 exp(-1.13eV/kT)C eq cr^sec" 1 . (24) 

1 i s 

Eqs. (21a) and (22) to (24) yield for the rate of converting uto Au g via 
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the kick-out mechanism to that via the Frank-Turnbull mechanism as 


r “ (dC g /dt) ko / (dC s /dt) FT = oexpC-S^leV/kT'/C®* 1 , (25) 

5 eq 

where a is a constant of the order of 10 . We do not kiow che quantity Cy but it 
can be seen from Eq . (25) t£at r decreases with a deriease of temperature unless 
the V formation enthalpy AH^ is equal to or larger than 3.71eV, which, judging 
from Eg. (21b), is unlikely. In Fig. 6 we show t* * calculated values of r(T)/ 
r(1000 C) assuming AH^ is in the reasonable range of l.b to 3eV. It is seen that 
the effetiveness of tne kick-out mechanism decreases much more rapidly with a 
decrease in temperature than that of the Frank-Turnbull mechanism. 

With this discussion, it becomes understandable vhy Au diffusion can appear 
as dominated by the kick-out mechanism (related to I) at higher temperatures, 
while by the Frank-Turnbull mechanism (related to V) at lower temperatures. It is 
expected that other metals may also exhibit a similar kinetical behavior in Si 
for experiments involving their diffusion; in- and out-diffusion experiments, 
precipitation and precipitate dissolution experiments. We suggest that by now 
the following rule should be kept in mind: in an experiment if it is found that 
I in Si appears as dominating the diffusion process of an impurity, it means that 
I exist as a species of Si point defect under thermal equilibrium conditions and 
the experimental conditions are suitable for I to appear as dominating, it does 
not mean that V do not also exist as a species of Si point defects under thermal 
equilibrium conditions; the reverse is also true. The diffusion experiments of 
Kitagawa et al. (41) cn Ni (dominated by V) and -f Mantovani et al. (42) on Ft 
(dominated by I) are good examples. 


4. ON THEORIES OF HIGH CONCENTRATION PHOSPHORUS DIFFUSION INTO SILICON 


Diffusion of jup-III and -V dopants shows several anomalous features that 
are quite thoroughl discussed in the literature (43-47). To summarize, 

Willoughby (45) and GUsele and Strunk (47) mentioned that the most prominent 
features are: (i) The emitter push effect, which is the extremely rapid diffusion 
of base dopant (B or Ga) of the doubly diffused npn transistor structure which 
results in an enhanced movement of the base-collector junction; (ii) The kink- 
tail structure of the in-diffusion P profile, which is resulted from the use of 
high concentration surface sources. This is shown in Fig. 7a with conveniently 
designated surface and tail regions. In the tail region the P diffusion rate is 
much higher than expected for isoconcentration studies or OED studies; (iii) 
Movement of hurried dopant layers can be extremely rapid if the surface P concen- 
tration is high. There exist ample evidences that P rich precipitates form during 
P in-diffusion (48-55), we therefore add to the above list by (iv) During P in- 
diffusion, the P concentration ir. the 31 surface region e ^ds its normal solu- 
bility limit. We consider that the following facts have be~n established: (i) 

The dopant diffusivity enhancement are caused by supersaturation of intrinsic 
point defects. The enhancement are much large** than that due to oxidations; (ii) 
The point defect supersaturations are related to high surface concentration of P, 
and, in some cases, of B or As (56-58); and (iii) Mechanisms involve disloca- 
tions in an essential manner, e.g., dislocation climb, are not the origin of the 
peine defect supersaturations. 
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Since untill recently it was widely believed that V are the dominant thermal 
equilibrium point defect species in Si, earlier models (59-61) put forth to exp- 
lain the P anomalous diffusion have been based on calculating the V supersatura- 
tions. Among these models, that of Fair and Tsai (61) gave a most complete treat- 
ment and was apparently able to calculate quantitatively the kink-tail structure 
in the P diffusion profile. Even if V were the only point defect species present 
in Si at high temperatures, the model still can not be regarded as satisfactory 
because of: (i) As pointed out by Kroger (62), the analytical development is que- 
stionable with respect to the relation between the total P concentration and the 
total carrier concentration on the one hand, and with respect to the derivation 
of the electrical field on the other; and (ii) As pointed out by Mathiot and 
Pfister (63), the V supersaturation was accounted for by a multiplication factor 
for deriving the fast diffusing tail region P diffusivity instead of solving the 
V diffusion equation. By correcting these problems Mathiot and Pfister (64) found 
that the model of F3ir and Tsai can not give a good qualitative description of 
the P diffusion profile. Here we mention two fundamental objections generally 
applicable to these models: ^i) Formation of P-rich precipitates (48-55) is not 
accounted for; and (ii) By now there also exist ample evidences that during P in- 
diffusion an I and not V supersaturation is present in the Si interior (55,64-67) 
but the V based models can only predict a V supersaturation. 

Recently, twc new models have been proposed. One is a vacancy percolation 
(modified E-center) model due to Mathiot and Pfister (68), which has also been 
extended to include an I-componen*; in the dif f usivities of both dopant diffusion 
and Si self-diffusion (69). The otner is a model due to Hu et al . (70) who pro- 
posed that a substantial part of P occupies interstitial positions which we shall 
call interstitial P model. Unfortunately, these two models seems are also not 
satisfactory. In the following (Sections 4.1 and 4.2) we discuss some essential 
points of these two models and give some detailed reasons for arriving at this 
conclusion. In section 4 . 3 we discuss a possible reason for P dift« ing into Si 
in a supersaturated concentration. 

4.1. The V Percolation Model 

By correcting the problems or the model of Fair ana Tsai and by incorpora- 
ting the idea of V percolation (71), Mathiot and Pfister were able to quantita- 
tively fit some P diffusion profiles (68) . Later on they were also able to extea 
the model to include an I-component in the dopant and Si self-diffusion coeffi- 
cients with an apparently even better success (69). Their earlier model (68) can 
be straight forwardly regarued as unsatisfactory since it is purely V based. \he 
incorporation of an I-conponent in the dopant and Si self-diffusion coefficients 
(t 9) is a goou step forward, but the model is not satisfactory because: 

(i) Formation of P-rich precipitates is not accounted for; 

(ii) They o concluded that G^-O.l to 0.3 for the dopants B, As and P from 900 to 
1200 C for low and high concentration diffusions. For the cas* of higher 
temperature experiments, we expect that local dynamical equilibrium between 
I and V to be reached, this can only result in a net V supei saturation in 
the Si interior. This is in contradiction with experimer.nl results (55,64- 
67) that can only be explained by a net I supersaturation; 

(iii) The finding that G^.<0.5 for P, which requires that the observed ORD oi P 
for (111) Si wafers oxidized for long times be explained -without assuming a 


266 



net V supersaturation (as Francis and Dobson (72) and the present autnors 
(8) have assumed), is in contradiction with the Sb OED ic&uii (10,19) which 
was observed under an almost identical experimental ooncition for observing 
the P ORD. A qualitatively consistent interpretation o these results is 
not obtained unless G >0.5 for P and G^<0.5 for Sb and under the experimen- 
tal conditions a net V supersaturation exists. 

(iv) The ooservation of ORD of B at 1200°C bv Hill (73) can not he explained by 
their finding G <0.5 for B, since in that experiment a 6% HC1 was used. By 
analysing the OSF data of Shiraki (21) the present authors nave shown that 
for •'his experimental condition a net 7 should have been injected (7,10). 

If Gj<0.5 f^r B, then it should have resulted in an OED instead of the 
observer OPD. 

The present authors believe that the V percolation idea is a very important 
point for which the full significance seems have not yet bee & realized. This sub- 
ject Js certainly a nice c ntribution of Mathiot and Pfister is weil worthy 
of some discussions. To invoke the V-mechanism for the diffusion of a substitu- 
tional impurity species which exhibit a higher diffusivity than that of self- 
diffusion of the host crystal atoms, it is necessary to assume that the diffusion 
proceeds in the form of i-V pairs (i denotes a substitutional impurity atom), and 
a pair has to partially dissociate during its rnigrat n process. In Si the sad- 
dle point for the i-V pair dissociation process corresponds to having a V on the 
third nearsi neighbour position of i. The remaining binding energy between i and 
V determines the cohesion of the pair. For the case of having a ver^ high impu- 
rity concentration, the impurity atoms effectively form a network cr cluster in 

the crystal lattice in the sense that now each impurity atom /ould have some 

other impurity atoms situated at its fifth or even closer neighbour positions. 

The i-V pair diffusion process changes now since upon reaching the saddle position 
the pair partial dissociation process, the V can also be at a third neighbour 

position of another i. Now V can migrate to the other i atom with a potentfcl 

barrier smaller than that needed for the isolated i-V pair case. This means with 
a sufficiently high impurity concentration, V can freely diffuse through the i 
network or cluster, which does not only enhance the V diffusion but also the i 
diffusion. This is the idea used by Mathiot and Pfister in their models (68,69). 
We now mention a few points associated with the percolation idea that have not 
yet been explicitly considered: 

(i) We believe that the percolation idea is equally applicable to an I-mecha- 
nism; 

(ii) In the presence of a high impurity concentration, and hence the percolation 
network, we believe it Is sensible to also suggest that the point detect 
formation energy is also lowered. This means the thermal equilibrium point 
defect concentrations are increased because of the presence of a u igh impu- 
rity concentration (not to be confused with that due tc electrical charge 
effect) ; 

(iii) The presence of the percolation network necessarily changes the i-I and i-V 
pairing potentials and hence G for low and high imparity concentration 

cl es are in princ. ^le of different values. 
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4.2. The Interstitial Phosphorus Model 

Hu et al. (70) have recently criticized nost of the V models (59-61) and 
proposed a new model to remedy the situation. An important contribution of their 
work is that formation of P-rich precipitate during P in-diffusion is very much 
emphasized and therfore in the model this is accounted for and hence removing 
one major sho; comings common to all other models. They have written down the 
basic equations but not yet performed the needed calculations and it is there- 
fore not possible to see whether a net I supersaturation in the Si interior can 
result from the model or not . An examination of the basic assumptions involved in 
this model has led the present authors to believe that the model in its present 
form is probably also not satisfactory. Their major assumptions are: (i) A por- 
tion of the P atoms in Si occupies the substitutional sites (P ) while the rest 
oc^ >^es interstitial sites (P^); (ii) Mutual conversion betweln P and P^ is 
Vtfi y .»iov; and (iii) Both P and P are shallow donors. Assumption 5 (i) is\ieeded 
since they invoked a two stream diffusion mechanism as the a priori reason for 
obtaining the kink-tail diffusion profile. The best way tc see what is a two 
stream diffusion is to imaging that there are two distinct chemical species A 
and B diffusing *nto Si simultaneously but otherwise independently and we are 
looking for C A+ g*C A +Cj,, see Fig. 7b. Now there is only one P species chemically 
but in Si they can occupy different sites, e.g., P and P.. In so far as diffu- 
sion is concerned, and P^ can behave independently. Assumption (ii) is needed 
since if P and P are converting into each other rapidly enough then their indi- 
vidual identity will be lost and a kink-tail structure in the P concentration 
profile can not be obtained. Assumption (iii) is needed since atomically disso- 
lved P atoms in Si is known to be electrically active. The assumption that a sub- 
stantial portion of the P atoms occupies inters tial sites is very different from 
that invoked conventionally: P atoms occupy substitutional sites only. By defini- 
tion, P diffuses primarily via an I-mechanism. Thus, there will always be a sub- 
stantial number of P atoms on interstitial sites whether they are migrating or 
not. In the converntional picture, P atoms can diffuse via an I-mechanism and 
hence those migrating P atoms are on interstitial sites temporarily but there 
will only be a negligibly small fraction of P atoms that ate migrating at any 
instant in time and hence for practical purposes all P atoms are occupying subs- 
titutional sites. 

In the following we present four arguments which jointly constitute as a 
strong indication that the above mentioned assumptions involved in the model of 
Hu et al. do not form a realistic enough basis for regarding the model as satis- 
factory (admittedly, though, each argument individually may be only constituting 
as a weak indication) : 

(i) The fact that there is but one shallow donor level associated with P doping 

is generally accepted. While it is legitimate to postulate that P is also 

a donor, it is difficult to believe that the donor is also a shollow one, 

and, in particular, its energy level needs to be almost identical to that 

due to P ; 

s 

(ii) High concentration P diffusion into Si can induce the formation of a dislo- 
cation network to accomodate the Si lattice parameter change in the diffused 
layer due to the incorporation of P atoms to a large concentration. Analyses 
of such dislocation networks indicated that the P diffused layer has a smal- 
ler lattice parameter (74) which is consistent with the knowledge that P 
atoms are all on the substitutional sites (the ratio of the covalent radii 
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of P atoms to Si atoms is 0.94)* In order to explain the P diffusion pro- 
files by the two stream mechanism, it needs to be assumed in some cases 
C can be almost equal to C p . If this is so, it seems that by now there 
should be some indication that the opposite, i.e., the P diffused layer has 
a lattice parameter larger than that of the undiffused Si, is true; 

(iii) The requirement that P^ and P are converting into each other very slowly 
is an extremely stringent condition. Considering the diffusion times invo- 
lved in a typical experiment, a simple estimate would show that the energy 
barrier against the P^-P conversion needs to be higher than 3eV and quite 
likely in the neighbourhood of 4 to 5eV. This is the same energy barrier 
prevent the interstitialcy pair from dissociation, and hence the pair needs 
to be regarded as very stable. If so, it is not easy to understand why * 
can diffuse so fast, which is also required in the model; 

(iv) In an Au gettering experiment, Lecrosnier et al. (75) found that the Au 
concentration in the P diffused layer developed a kink-tail profile similar 
to the P profile. This phenomenon is consistent with the interpretation 
that all P atoms occupy substitutional sites but not equally consistent 
with the interpretation that a substantial portion of P atoms occupies 
interstitial 3ites. The interstitial P model requires that the surface re- 
gion P atoms be mostly P (diffusion slow) while the tail region P atoms be 
mostly P (diffusion fasl). We first mention that in the Au gettering expe- 
riment all the gettered Au atoms should be Au since normally C e< *>>C^ holds 
and since if they were Au^ then we expect tfie P and Au concentration pro- 
files follow a complementary rather than a similar distribution pattern 
since the charge states of the P and Au atoms should result in a repelling 
force (P is a donor in spite of whether in substitutional or interstitial 
positions, and, Au is also suppose to be a donor, i.e., Au; . Assuming 
that Au atoms are indeed Au atoms, then the best interpretation is that P 
atoms are all P atoms s ^ce now charge compensation occurs between P + and 
Au^ (presumably Au ) atoms and atomic size compensation also occurs (cova- 
lent radius of Si atom is larger than that of P but smaller than that of 
Au ) . If the P atoms are then we expect a large compressive stress to 
develop In the Si which tends to result in a complementary distribution 
pattern between the P^ and Au concentrations. The worst case is that both 
P and Au atoms occupy interstitial sites since then it will result in a 
compressive strain as w.ll as a charge repelling situation. 

4.3. Physical Reasons for ? Diffusing into Si in Supersaturation 

While the kink-tail diffusion profile of P has been recognized as an out- 
standing feature for a long time, it is only recently that the equally outstan- 
ding (and undoubtly related) feature of P diffusion into Si in supersaturation 
is also recognized. Why is it so? Hu et al. (70) provided one possible reason 
which is dynamical in nature: P diffusion (using a P-rich oxide source) is usua- 
lly carried out in a partial oxygen ambient and hence oxidation, i.e., a chemi- 
cal reaction, occurs. They suggested that the excess chemical potential respon- 
sible for P entering Si exceaaxcg its solubility limit derives from this chemi- 
cal reaction (oxidation). However, this may not be the primary reason because 


*This problem was brought to the attention of one of the \ ^esent authors (TYT) 
by Dr. T. E. Seidel. 
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of: (i) The I supersaturation associated with high concentration P diffusion far 
exceeds that due to an oxidation of Si, and, in particular, this seems to be 
true at very low temperatures, e.g., 750 C, for which an oxidation of Si ceases 
to function; (ii) There are some indications that P can diffuse into Si in a 
supersaturated amount for cases that the annealing ambient is inert. Nobili et 
al. (54) concluded that the electricallg-actiYe P concentration limit (aresuma- 
bly the solubility limit*) is about 4x10 cm at llOO^J, bu^ at the sam^ tempe- 
rature P can enter Si to a ^roentration of about 1x10 cm J in an dr. T vt=- in 
annealing using a P-rich oxide source (76). 


As an alternative possibility, we suggest a mechanism due to phase-equili- 
brium requirements to explain the phenomenon of P diffusing into Si in supersa- 
turation. This mechanism is readily understandable if the materials involved 
constitute as a binary alloy system. Therefore, for our present purpose we 
discuss the case that the source of diffusion is gaseous P (P ) at the standard 
presure instead of the more commonly used P rich oxide. In Fi§. 8 we show the 
phase diagram of the Si-P system and an appropriate free energy diagram in the 
temperature range of our interest. The significance of Fig. 8 is that, under 
phase-equilibrium conditions: (i) the two primary phases, a (Si containing ato- 
mically dissolved P) and y (P ) , and the binary intermediate phase SiP, can each 
exists alone; (ii) The a phasl and SiP can coexist; (iii) SiP and P can coexist. 
This means for a diffusion experiment using P as source, to start §ith, the 
materials (Si and P ) do not exist in a phase^equilibrium state. Hence, the ma- 
terials will proceel to reach the appropriate phase-equilibrium states by crea- 
ting coexisting phases permitted by the phase-equilibrium conditions, i.e., the 
above discussed possiblities (ii) and (iii). We coined the term that, before the 
phase-equilibrium state is reached, the materials are in a diffusion-transient 
state. A primary phenomenon associated with this transient state is that P will 
diffuse into Si in excess of its normal solubility limit which is defined per 
coexistence of the a and SiP phases. Refer to Fig. 8b, the normal P solubility 
in Si is represented by concentration A. The material near the Si surface is now, 
however, not coexisting with SiP but with P which gives a P solubility in Si 
represented by concentration B. Since B>A his to hold, we see P enters Si at a 
concentration exceeding its solubility limit A. That is, by the very fact that 
now Pg and Si coexist, the material Pg directly exerts a chemical potential on 
Si to maintain a P concentration in excess of the P solubility deverived per co- 
existence of Si and SiP. If kinetical conditions permit, SiP precipitates can now 
form inside the Si crystal or alternatively a continuous layer of Si? can form 
at the Si surface. Before a sufficiently thick continuous SiP layer forms, the 
supersaturation phenomenon will persist. 

Per this discussion, we see that there is a thermodynamic reason for P 
to diffuse into Si in supersaturation. For experiments for which a P rich oxide 
layer is used as source, the argument is also applicable. The analysis leading 
to this conclusion is, however, considerably more complicated than the binary 
example just discussed. This complexity arises because a Si-0-P ternary situation 
is involved. We have carried out the analysis for this ternary situation and the 
details will be given elsewhere (77). In that analysis we have also shown that 
the 0 solubility limit is also changed. There exists one experimental result (76) 
indicating that this seems to be the case. 
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5. CONCLUDING REMARK 


In this paper some essential points associated with oxidation effect - : and 
Au diffusion in Si are discussed. Progresses in these areas in the last few years 
seems to have established beyond reasonable doubts that I and V coexist in Si 
under thermal equilibrium conditions. Within the frame work of coexisting I and 
V species, some new understanding is now emerging. We look forward for more pro- 
gresses in our understanding of the nature of point defects and their roles in 
the materials transport properties in Si in the near future. 
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Fig. 1. Prediction of dopant dif- 
fusivity changes due to an oxida- 
tion. The model assumes that I and 
V coexist in Si in thermal equili- 
brium and in oxi ition at high 
temperatures ana that during the 
oxidation I and V attained local 
equilibrium. 
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Fig 2. Quantitative fitting of aval 
lable Sb ORD data from (100) wafers 
to the calculated Gj s 0.02 curve per 
Eq. (8). Unfilled circle: Mizuo and 
Riguchi [18]; filled circle: Tan and 
Ginsberg [19]. 




10 A /T 

Fig. 3. The recombination bar- 
rier limited time constant, 

T bar» * or 1 and V to reacl * 
dynamical equilibrium as a 

function of reciprocal tempe- 
rature calculated per Eq. (ll). 


Fig. 4. Gold concentration 
profile after a lhr diffusion 
into dislocation-free Si at 
900°C [14]. Dash line: an 
erf c-function fitting per Eq . 
( 19 ) ; solid line: fitting 

per Eq. (20). 
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Fig. 5. DjC®** and vs. recip- 

rocal temperature. 



Fig. 6. The function r(T)/r(10C0 C) 

. vs. reciprocal temperature, r is the 
ratio of the efficiency of converting 
Au into Au due to the kick-out me- 
chanism to that due to the Frank- 
Turabull mechanism. Numbers above 
each curve are the possible AH^ 
values . 



Fig. 7. (a) Schematic dra 'ng of the kink-tail structure of P profile, 
which, hypothetically, may be obtained by adding (b) profiles of phy- 
sically distinguishable A and B atoms diffusing independently. 
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Si A B p 

(b) 


Fig. 6. The Si-P binary system phase diagram (a) and an appropriate free 
energy diagram (b) showing the thermodynamic origin for P to diffuse into 
Si in excess of its normal solubility limit. The solubility limits A and 
B are obtained from the usual common tangent constructions per coexistence 
of the o phase with the SiP and Pg phases respectively. 



DISCUSSION 


BOSK: On one of your viewgraphs yon showed during the phosphorus diffusion, 

something like 40% of it is electrically active. I wonder if you might 
coament further on what physical form does it take? 

TAM: It takes precipitated form. A small part of it, about 15%, is identified 
as SiP structure. About 85% of them are only observed in terms of no more 
than 15-A-sised precipitates. They are too small to be identified by 
transmission electron microscope. Therefore they are only inferred to be 
Sip precipitate, but not proven. 
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